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Abstract 



Continuing our effort to build a consistent power counting for chiral nuclear effective field theory, we discuss 
the subleading contact interactions, or counterterms, in the singlet channels of nucleon-nucleon scattering, 
with renormalization group invariance as the constraint. We argue that the rather large cutoff error of 
the leading amplitude requires 0{Q) of the low-energy approximation to be no n- vanishing, contrary to 
Weinberg's original power counting. This, together with the ultraviolet divergences of two pion exchanges in 
distorted- wave expansion, leads to enhancement of the ^5*0 counterterms and results in a pionless theory- like 
power counting for the singlet channels. 
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I. INTRODUCTION 



The power counting is one of the essential ingredients of any effective field theory (EFT), which 
not only keeps track of an infinite number of operators and Feynman diagrams but also estimates a 
priori the neglected contributions for a given order. Naive dimensional analysis (NDA), a corner- 
stone of Weinberg's original power counting [l| (WPC) for few-nucleon systems, is often employed 
to assess the size of coupling constants: each derivative on the Lagrangian terms is always sup- 
pressed by the underlying scale of chiral EFT, Mhi. Though phenomenologically successful 
WPC has been shown to be inconsistent with the principle of renormalization group (RG) invari- 
ance 7-l3], especially in the triplet channels where the singular attraction of one pion exchange 



(OPE) calls for modifications to WPC at as early as leading order (LO) [lO|. The issues of RG 
invariance are less acute in the singlet channels since the LO amplitudes of WPC for these channels 
are indeed RG invariant 0, lol. [isl. [l^ . But a modification to WPC for the subleading countert- 
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to be necessary. Following our 



erms in the singlet channels has been argued in Refs. 

investigation of the triplet channels [o] > we use RG invariance as the guideline to study the 
subleading counterterms of the singlet S and P waves. 



The EFT expansion of the T-matrix at low energies has the quintessential form, 

where Q denotes generically external momenta, n the counting index, Mio low-energy mass scales 
and Fn{x) the non-analytic functions from loop integrals. Even though the nonperturbative uni- 
tarity requires any nonrelativistic, nonperturbative T-matrix to scale as Q~^, we choose to label 
LO as 0(1) so that one does not need to change the standard chiral power counting for irreducible 
pion exchange diagrams, e.g., OPE is 0{1). Subleading orders of the EFT expansion are labeled by 
their relative correction to LO, i.e., next-to-leading order (NLO) by 0{Q/M\y\) or 0{Q) for short, 
and next-to-next-to-leading order (NNLO) by 0{Q'^ jM^-^ or 0{Q'^), and so on. 

In any EFT calculations employing a ultraviolet (UV) momentum cutoff A, the cutoff 
independence — RG invariance — of the T-matrix is usually imperfect at a given order. There 
usually exists a residual cutoff dependence, or a cutoff error, of the T-matrix with a large but finite 
cutoff, which vanishes though as A oo. Seen in this light, RG invariance needs a little more 
careful interpretation. Because (i) the cutoff error is part of the theoretical uncertainty at a given 
order and (ii) the theoretical uncertainty is, by definition, of the same order as the next-order EFT 
correction, the cutoff error must be smaller than or of the same order as the next-order correction. 
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This seemingly trivial statement constrains power counting in a non-trivial way when certain order 
is considered vanishing. Of our interest is WPC, in which the 0{Q) corrections have long been 
deemed to be zero; thus, the theoretical uncertainty for LO is considered by WPC to be 0{Q'^). 
It follows that the LO cutoff error should vanish at least as fast as Q^/A^: 

(g; A) -rW(g;oo) (2) 



While this is the case for the triplet channels [3, IS 



21[, we will show that it is not for ^Sq, which 



forces us to modify WPC for ^Sq at subleading orders even though its LO satisfies RG invariance. 

If a counterterm is not required by RG invariance at 0{Q^) but is counted in NDA, we 

will follow NDA to power count that counterterm. In other words, we do not minimize the number 
of counterterms at a given order using RG invariance as the criterion 



13, 



22l |. The rationale for 



this is as follows. RG analysis does not study the degrees of freedom that are not built in the 
effective Lagrangian, such as non-Goldstone bosons and/or very heavy excited baryon states. To 
avoid underestimating the contributions of these degrees of freedom to counterterms, we set the 
minimal size of counterterms as the one given by NDA. 

We will treat subleading potentials as perturbations on top of the LO T-matrix, which is 
nonperturbative iteration of the LO potential. This is sometimes casually called perturbative 
renormalization. Underlying this approach is the point of view that power counting should be 
done at the level of physical observables, in our case, the on-shell scattering amplitude. Infrared 
enhancement due to nucleon intermediate states should be incorporated into the power counting [23l . 

rather than be used as the pretext to settle for WPC. However, this is not to say that a consistent 
power counting for the nonperturbative approach, in which the full iteration of the whole potential 
is performed, cannot be found. For the development in this direction, we refer the reader to 



Refs. 



Our paper is structured as follows. In Sec. [Ill the subleading counterterms are classified into 
three categories according to the diagrams that drive their RG evolution. After a short review 
of LO in the singlet channels, we discuss in Sec. IIIII the rather large cutoff error of the LO ^Sq 
amplitude and the consequence of that for power counting. We analyze in Sec. IIVI how the non- 
vanishing 0{Q) counterterm and two pion exchange contribute to 0{Q'^). This is followed by a 
discussion and a conclusion offered in Sec. |Vl 
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II. EVOLUTION OF SUBLEADING COUNTERTERMS 



The LO T-matrix, T^^\ arises from the full iteration of OPE (in low partial waves) and a set 
of counterterms that ensure the RG invariance, which means that T^^\k;A) is independent of 
A when A ^ A;, where k is the magnitude of the center-of-mass momentum. Subleading orders 
are given by perturbative insertions of higher-derivative counterterms and/or irreducible multiple 
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28l | is far more intricate, we expect 



pion exchanges into LO. Although renormalization at LO 
to have a better visualization of renormalization at subleading orders by forming a fairly simple 
correspondence between a loop diagram and the counterterm to subtract its UV divergences, much 
like ones that exist in perturbative EFTs. In Wilson's language of RG analysis, this counterterm 
is the one that "evolves" most significantly when the cutoff of the loop diagram is rescaled from A 
to a smaller value. A', but remains large in the sense A' ^ k. It will help our discussion to classify 
loop diagrams and their corresponding counterterms into the following three categories. 



A. Residual counterterms 

The first class of diagrams are actually those of T^^\ Of course the LO contact operators, by 
definition, are the counterterms to renormalize T^^^; however, in order to systematically remove 
the residual cutoff dependence of T^'^\ one must take account of the contact operators with more 
derivatives than the leading one. We call those higher-derivative operators the residual countert- 
erms for r^''). Even though they are not as important as the LO counterterms or OPE, they might 
be more important than multiple pion exchanges which start to contribute at 0{Q'^). Its important 
for us to find a way to estimate their sizes before calculations are carried out. 

The authors of Refs. [111. l2d. 1291] have attempted to analyze both LO and residual counterterms 
using the Wilson RG equation, with OPE as the only long-range force. This is a very difficult task 



and several assumptions were made in Refs. 
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2Q|, 



29|. Energy and momentum dependences of 



contact operators were assumed to be independent of each other. But we know that, when treated 
as perturbation on top of LO, they can be related by the equation of motion. The RG invariance 
of the off-shell T-matrix was imposed, though only the on-shell quantities need to be RG invariant. 
While this excessive requirement cannot be deemed wrong, one may be concerned that the resulting 
power counting demands more counterterms than necessary, only to ensure the RG invariance of 
the off-shell part of the T-matrix. What is most debatable is perhaps the existence of a fixed-point 
solution to the RG equation, around which the power counting is obtained. Although this appears 
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to be reasonable in the singlet channels, it is clearly at odds with the running of counterterms in 



the attractive triplet channels, where a limit cycle-like behavior was observed 

Instead of performing a comprehensive analysis of the RG equation, we impose a narrower 
definition for residual counterterms: the correction brought by the residual counterterm is of the 
same size as the cutoff error: T^'^Vfc; A) — T^*^^ {k; oo). For example, the LO cutoff error in — ^Di 
is found to be 0{k'^ M^J'^ A'^^"^) [8, 9, 21]. Thus, the correction due to the ^Si residual counterterm 
is rated as ©(Q^Afj^^A-^/^)^ jggg ^i^g^^ ^j^g leading two pion exchange (TPEO), which is ©(Q^/M^j). 

It should now become apparent that we use residual counterterm merely as a mnemonic device 
to reflect the order of magnitude of the LO cutoff error. With this notion we can interpret the 
inequality ([2|) as follows: WPC requires the residual counterterm to be no more important than 
TPEO. It is indeed true for the triplet channels, but, as we will show, not so for ^5*0. 



B. Primordial counterterms 



Better known are the second class of diagrams: irreducible multiple pion exchanges evaluated 



171 ] are the contact operators 



in the plane, or free spherical, wave basis. Primordial counterterms 
necessary to subtract the divergences of these pion exchange diagrams. 

On the basis of "naturalness" , the power counting of the primordial counterterm should be the 
same as that of the pion exchange diagram in question, which is reliably handled by WPC. For 
example, TPEO is 0{Q'^) and its primordial counterterm, a second-order polynomial in momenta, 
is counted 0{Q^) as well. 



C. Distorted-wave counterterms 



The third class of diagrams are insertions of irreducible multiple pion exchanges into LO, that 
is, pion exchanges sandwiched between distorted waves — the LO wave functions, Vfc- The coun- 
terterms to absorb the divergences of these diagrams are called distorted-wave counterterms. 

Coordinate space provides the best stage for qualitative discussion of distorted-wave UV diver- 
gences. Any reasonable UV regulator will roughly separate radial coordinate into two parts: the 
inside, < r < A~^, and the outside, A~^ ^ < oo. Details of the regulator decide how sharp the 
separation is. When contact interactions are present at LO, such as in both S waves, the inside 
and outside part of the LO wave function, i/;™ and "0^"*) are respectively subject to the LO contact 
potential and OPE. The LO contact potential and OPE have different short-distance structures, 
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which result in different short-distance behaviors of i/'™ and V'fc"*- For instance, as shown in, e.g., 
Refs. jisl. [30I. tp'^^^{r) in ^Sq has an irregular component diverging like ~ 1/r near r ~ A~^, which 
differs very much from a free spherical wave. However, since the LO contact potential is always 
well defined upon regularization, the inside wave function -0^'^ is not drastically different from a 
free spherical wave. 

The distorted- wave matrix element of a subleading contact potential, (V^fcl^J^^l^fc)? is domi- 
nated by the integration over the inside region. Since V'aTI^) behaves similarly to a free wave at 
short distance, {'4'k\Vs'^^\'4'k) is expected to be as UV singular as its free-wave counterpart, i.e., 
(plane wave |y|"''| plane wave). 

On the other hand, the distorted-wave matrix element of a subleading long-range potential, 
(^fc|^"^|^fc)) is mostly decided by the integration of ip'^^^{r). The irregular component of 
if present, could make (V'fcl^^^lV'fc) more divergent than its free-wave counterpart. Therefore, the 
primordial counterterm that renormalizes a given multiple pion exchange in the plane wave basis 
may no longer renormalize the same pion exchange between distorted waves. If it does not, the 
distorted-wave counterterm will have to be more singular than its primordial counterpart, that is, 
will have more derivatives. This indeed happens to a toy model considered in Ref. 

16] and ^Sq [l^, among possibly other channels of A^A^ scattering. We will 
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3ll | and in 



reproduce in Sec. IIVI the distorted- wave enhancement of the subleading counterterms, which 
was first shown in Ref. \\\ . 

The distorted- wave multiple pion exchanges, (■0fc|V^"^|'0fc)) ^-re power counted as the same as 
the free-wave matrix elements because it has been established at LO that any number of insertions 
of yC^) — the LO potential — does not enhance or diminish the amplitude. Based on, again, 
naturalness, it follows that the distorted-wave counterterm is power counted as the same as its 
primordial counterpart, even though one of them may have more derivatives. 

III. LO AND ITS RESIDUAL COUNTERTERMS 

We consider first ^Sq. The LO amplitude is constructed by resumming which is a constant 
counterterm Cisq plus OPE, 

y(o)(g) = K(g)+C.5o, (3) 

where 

^-(9) = rl^x-3' (4) 
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with p' (p) the outgoing (incoming) momentum in the center-of-mass frame, q = p'—p,gA = 1-26 
and /,r = 92.4 MeV. We can redefine Cig^^ such that the point-hke piece embedded in OPE is sepa- 
rated from the Yukawa potential and rewrite long and short-range potentials for ^^o, respectively, 
as 

vy = -^^^, = c(°), (5) 

rriM q + 

where a^^ = g\mj\f/16-Kf^ ~ (290MeV)~^ and A-K/mj\f is a common factor of nonrelativistic Feyn- 
man amplitudes. Here we have dropped the subscript ^Sq to simplify the notation. C has been 
expanded in anticipation that the running of C(A) with respect to A could be modified at each 
order, 

C(A) = ^(0) (A) + C(^) (A) + • • • , (6) 

though the number of physical inputs to determine C remains one (or stated differently, the bound- 
ary condition for the RG fiow of C (A) remains fixed) . There is little dispute about the power 



counting of (renormalized) C 



24 



C.~^T^, (7) 



where Q ~ Mjo ~ a^"^. The relatively large size of a^^, compared with ttt-tt, is crucial for the 



singlet-channel success of KSW scheme 



23|, 



331 ] ■ in which OPE is treated perturbatively. But here 
we take the view point that is an infrared mass scale because (i) it is still smaller than Mhi 
and (ii) the perturbative Yukawa works less well when takes a larger value but still stays within 
the validity of chiral EFT. 

To study the UV divergences, we find useful the elegant machinery developed in Ref. Q]- First, 
we define the resummed Yukawa amplitude: 

Ty(/,p ; A;) = Vy{p',p) + t jP^Vy{p', f; k) J^^l^^^ , (8) 

where E = k'^ /rriN is the center-of-mass energy. The superscript "A" only serves as a reminder 
that every integral is in principle regularized. While the LO P and higher wave amplitudes are 
given solely by the resummed Yukawa, the LO 5'-wave amplitude needs summing up insertions of 
C^^^ to all orders, which is eventually given by Q] 

T"- '{p ,p;k) =Ty{p ,p;k)+ i — , (9) 

(Cy^i) -Ik 
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x(p; k) = 
Ik = 



FIG. 1: Diagrammatic representation of xiP'i ^) ^'^d Ik- Here the solid (dashed) hne represents the nucleon 
(pion) propagator, and the crossed circle docs not represent any interaction. 

where 

X{p; k) = l+ / j^— , (10) 

Figure [T] shows the diagrams that, when resummed, represent and x(pj k). The power counting 
of Ik and x(^i fc) will follow, e.g., that of the first diagram of their resummation series: Ik ~ '"^^ 
and x{k'i k) ~ 1- 

Although we cannot calculate analytically Ik, xip'^k) or Ty, the dominant UV divergences can 
be captured by noticing that every insertion of Vy — combined with the Schrodinger propagator — 
suppresses the UV divergences by 1/A Q]- Therefore, the most significant cutoff dependences of 
Ik are in the first two diagrams of the lower row in Fig. [TJ the first is linear and the second is 
logarithmic in A, 



—Ik = — {lA + I^)+m,-^,rnl)^ + o(^) , (12) 



niN 
with 

47r / A \ 

/A = /3oA + /3iA^^ln — . (13) 

Here I^ is the finite part, k^^ = aT^ni^ and j3i are dimensionless and depend on the details of the 
UV regulator. While /5o and /3i are numbers coming out of the first two diagrams, the 1/A cutoff 
dependence, /32, receives contributions from all the diagrams; therefore, (32 is a nonperturbative 
function of k, k^^ and m^. In the spirit of keeping track of 1/A cutoff dependences, one finds that 

x{k;k) = x§ + 72ik,^.,ml)^ + o(^^^ , (14) 

where Xk finite part and 72 is another dimensionless function. As for Ty, we notice that its 

residual cutoff is only 0(K7rfc^/A'^). 



We would like to have available the coordinate space form of the LO wave function, which will 
be useful in analyzing the distorted- wave counterterms for TPEs. Since they are somewhat out of 
the main line of our discussion, the relevant expressions concerning the wave function are relegated 
to Appendix lAl 

Introducing the renormalized coupling, Cr, such that 



we can rewrite the on-shell T-matrix as 



C^1 = (CW) '-/a, (15) 



r(0) {k, k- k) = Ty (k, k-k) + — . (16) 



However, one can go 



a. 



As A — )■ oo, this is no more than reproducing one of the results in Ref. 
further and infer from the cutoff error some information about the subleading ^Sq counterterms. 
Using Eqs. ([T2]) and ([H]), power counting d?]), Ik ~ "^^J^ and x{k',k) ~ 1, one finds, as promised, 
that the cutoff error is ©(fc^/AfioA): 



i?2 



„ 72 k"^ 

R -r + P2- 



(17) 



Therefore, the theoretical uncertainty of T^'^'f must be 0{Q), or equivalently, the residual countert- 
erm, the two-derivative ^5o contact operator D /2{p''^ +%P'), is 0{Q). On the other hand, since pion 
corrections do not start until ©(Q^), i.e., TPEO, a non-vanishing 0{Q) can only be one insertion 
of the D term: 



where Vg'^ denotes the short-range part of the 0{Q) potential. 



's 
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(%|yWli5o) = cW + (19) 



and G is the Schrodinger propagator. 

If the loop integrals in x(^; k) (fTOl) and Ik (fTTI) are dimensionally regularized, the LO cutoff 
error (jl7p vanishes. But it is a model-dependent thinking to count on a specific regulator to be 
superior and to ignore the cutoff error that may arise with other regulators. 

Higher partial wave amplitudes are decided by Ty alone. The quite small cutoff error of Ty, 
C'(K7r^^/A^), means that the residual counterterm for ^Pi is no more important than 0{Q^) and 
that WPC does not need to change for ^Pi: 

{^Pi\V^^\^Pi) = 0. (20) 



Through the steps shown in Appendix [Bl one can rewrite the 0{Q) ^Sq amphtude (jlSp in a 
more comprehensible form: 



(21) 



(C(o))2 [(C(o))- 

where V^^\x) is V^^^ represented in coordinate space. Since C^^^ does not incorporate any new 
physical information, it is at our disposal to choose the value of C^^^ so long as it helps renormal- 
ization. With 

D(0) 



we choose the value of C^^^ such that 

cm 



:Cl-mNDRV^''\0) = 0. 



(22) 



(C(0))2' 

Now we can express in terms of renormalized quantities: 



(23) 



^(1) _ Dr Pxf 



(24) 



It is obvious that, for to be 0{Q), the scaling of Dr must be 

47r 1 

Dr , (25) 

in comparison with Cr ([7]). The interpretation of becomes particularly simple at the chiral 
limit where Ty vanishes: the D term plays the role of the effective range. 



IV. 0(g2) 



S Wave 



At 0{Q'^) there are two insertions of and one insertion of each of TPEO (denoted by V2^^) 

(2) (2) 

and Vg , where Vg , before any higher-derivative counterterm is considered, includes at least the 
0(Q2) corrections to C(A) and D{k): C(2)(A) and L>(i)(A). 

Two insertions of Vg'^ include integrals involving the LO interacting Green function: 
G{1 + T'^^^G). The contributions of ^ and two vi^^'s are eventually summed up as 



(0)^ 



"k = 



if2 



p(2) ^ 



+ 



Dl C7(o) 



Cl {Cr' - /f ) Cr Cr 



(26) 
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where 



A = —mi\fD 




(0) 



ruN 



y(o)(o) 



5f (0) 



+ 



f<2 



2 ' 



^2 



(27) 
(28) 



with the regularized Dirac function 6)^ (f) defined by Eq. (IBOh . For the related computational 
details of the above equation (and Eqs. ([29]) and ([30l) ). we refer the reader to Appendix|Bl The first 
term in the bracket of Eq. (j26p does not bring more information than T^^) (j24p : it merely restores 
the unitarity up to 0{Q'^). The second term has new structure, which becomes more apparent at 
the chiral limit where, since Vy vanishes, it resembles the shape parameter of a contact-only theory. 
With power counting d?]) and (j25p and Eqs. (jl3p and (jlSp . one sees that the "shape parameter" 
term is C'(Q^/M^;(Mio + /3oA)). We cannot relegate it to 0{Q^) even though it vanishes at A — ;> oo 
because (3q is not necessarily 0{1) or larger. In fact, if dimensional regularization is used /3o is 
zero, which makes the shape parameter term 0{Q^) [^l- To control this piece of regularization 
dependence, we must add a four-derivative counterterm E'^^^p''^p^ to v!?^ . One insertion of the E 



term yields 



i?2 



iriN 



y(0)(0) 



(29) 



Before considering TPEO, we remark that the other four-derivative term S'{p'^ when 
treated as perturbation, is redundant: 



.(2) 



A;2 -m^y(°)(0) 



mivVVW(0)| 



(30) 



With Eqs. (I24p and ()29p . T\ can be expressed as a combination of the C and operators. A more 



general argument is of course the field redefinition inspired by the nucleon equation of motion [28| . 

The analytic part of {q) is a second-order momentum polynomial, i.e., its primordial coun- 
terterm, which, when projected onto ^Sq, is nothing but the D term. The insertion of the D term 
into T^'^) is shown in Sec. IIIII We now consider the matrix element of the non-analytic part, which 
diverges as r — ?• in coordinate space: {r) ~ l/(M^jr^). This was first shown in Ref. 1^ . 
though in a slightly different notation than ours. 



This is perhaps most readily done in coordinate space where V^^^^ is diagonal 



n(2) 



47r / drr'i^l{r)vl\r), 



(31) 
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where ipk{f) is the LO "^5o wave function. The "outside" part of ipk{i") {^~^ ^ f) is subject to the 
Yukawa and is dominated by the irregular solution Hkir) ()A3p at short distance, which diverges 
as 1/r near r ~ A~^. As a consequence, the UV divergence of Tg^'' is captured by the integration 
of the outside wave function from any infrared length down to r ~ A~^, 



of Kjt/A and have at most logarithmic dependence on A. While poA"^ and pik'^h? can be respectively 
subtracted by C^^) ([27]) and D^^") ([28]), the divergence proportional to p2k^\iiA needs E^^^i ([29]) 
to cancel. That is, the ^5*0 distorted- wave counterterm of TPEO — the E term — has two more 
derivatives than the primordial counterterm — the D term. 

Using the fact that a multiple pion exchange with 0{Q"' /M^-^) correction to TPEO behaves as 
l^^5+n gjjQj.^ distance and repeating the above procedure, we can eventually conclude that, for 
any multiple pion exchange, the ^Sq distorted- wave counterterm is a momentum polynomial with 
two more powers than their primordial counterpart. 

For the sake of consistency, we have kept track of residual cutoff dependence at each order and 
have assured ourselves that they are not lower-order than next-order correction. Besides being 
an increasingly cumbersome task, it is less important to show explicitly the cutoff error at each 
order; the power counting is much less demanding on the RG invariance because no certain order 
is considered vanishing beyond 0{Q) thanks to multiple pion exchanges. 

We have seen two motivations to promote the E counterterm to 0{Q'^): (i) to control the 
regularization dependence of two insertions of the D term and (ii) to absorb the distorted-wave 
UV divergences of two pion exchanges. 

Unfortunately, the integral in Eq. (|32p cannot be evaluated analytically; thus we cannot express 
the full 0{Q'^) amplitude in terms of the previously defined renormalized building blocks. But the 
structure of Vs at 0{Q'^) will suffice in the numerical calculations carried out later: 




(32) 



where Af/C^^^ is RG invariant (see Appendix|A]) and "F.T." refers to finite terms, pi are functions 



(33) 



B. P wave 



The distorted- wave counterterm for TPEO in Pi is the same as the primordial counterterm be- 
cause, without an irregular component, the LO P-outside wave solution cannot make the distorted- 
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FIG. 2: (Color online) ^5*0 phase shifts as a function of laboratory energy. The red dots are from the 



Nijmegcn PWA [Sj]. The dark green (hght blue) band is the LO {0{Q)) EFT resuh with A = 0.5 - 2 GeV. 
The dashed, dotted and dot-dashed lines are 0{Q^) with A = 0.5, 1 and 2 GeV, respeetively. 

wave counterterm more singular than the primordial one. It follows from this, combined with the 
observation that the residual counterterm for ^Pi is not larger than TPEO, that WPC does not 
need to change for ^Pi; a single P-wave counterterm is what is needed for 0{Q^) and 0{Q^), 



(ipi|yj2,3)|ip^ 



(-(0,1) „/„ 
^ip-^ PP- 



(34) 



Numerics 



We compare our EFT calculations with the Nijmegen partial wave analysis (PWA) j34l ]. The 
expressions for the delta-less TPEs from Ref. are adopted here. Plotted in Fig. [2] are ^5*0 phase 
shifts versus laboratory energy, Tiab- The LO curve is fitted to the PWA at Tiab = 5 MeV, and the 
PWA points at Tiab = 25 and 50 MeV are added to determine D and E, respectively, at 0{Q) and 
O(Q^). A good reproduction of the PWA is achieved up to Tiab — 100 MeV, which translates into 
k ~ 200 MeV. The sharp momentum cutoff is used in the calculation. 

Unlike in the triplet channels, the analytical arguments for renormalizability in the singlet 
channels are quite solid. So it is less crucial to examine numerically the cutoff (in) dependence of 
the EFT amplitudes. However, it is still reassuring to see that the 0{Q'^) curve with A = 1 GeV 
is closer to the A = 2 GeV curve, suggesting the cutoff independence for large A. 

Although WPC is intact for ^Pi, we plot in Fig. [3] ^Pi phase shifts, for completeness. The 
cutoff independence is rather trivial for ^Pi; therefore, only A = 1.5 GeV is employed. Since 
going to 0{Q^) in ^Pi is much easier than in ^5*0, we include 0{Q^) results as well, in which the 
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FIG. 3: (Color online) ^Pi phase shifts as a function of laboratory energy. The red dots are from the 
PWA Q. The dashed and dotted lines are respectively the LO and 0{Q^) EFT curves. The dot-dashed 
(set I) and solid (set II) lines are 0{Q^) curves with different sets of irirNN seagull couplings (see the text 
for explanation). 



subleading TPE (TPEl) contributes. There is only one counterterm up to ©(Q^), Cip^ dS]), which 
we determine by fitting to the PWA point at Tiab = 50 MeV. 

TPEl has crucial dependences on the ttttNN seagull couplings, Cj, that have chiral index v = 
1. We show 0{Q^) EFT curves, respectively, with two commonly-used sets of Cj (in unit of 



GeV"^): (I) the dot-dashed line with ci = -0.81, 
solid line with ci = —0.81, C3 = —3.4 and C4 = 3.4 [6 



4.7 and C4 = 3.4 [351] and (II) the 
I, [1^. The impact of the uncertainties of 



Ci is significant beyond approximately 50 MeV. Since the uncertainties of q have their roots in 
slow convergence of the delta-less description of vrA^ scattering, we expect that the delta-ful pion 



exchanges |2|. l36l. l37l|. with the irNA low-energy constants determined by the vrA^ scattering data 
around the delta peak will improve the convergence of chiral NN EFT. In fact, aside from 

the open issues of power counting counterterms, the delta-ful nuclear forces have been shown to 
achieve a more ra pid convergence in the two-nucleon Q, [4^ and, on a more qualitative level, in 
the three- nucleon [43| sectors. 



V. DISCUSSION AND CONCLUSION 



We have studied how RG invariance constrains in the singlet channels the structure of sublead- 
ing counterterms, with S and P waves as the examples. Our analysis shows a hierarchy of ^Sq 
counterterms that resembles the pionless theory while WPC remains unchanged in ^Pi and higher 
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singlet partial waves. 

To facilitate the discussion, the subleading counterterms are classified into three categories ac- 
cording to the loop diagrams that drive their evolution (see Sec. |Tl]). The residual counterterms 
eliminate the small cutofi' dependence of the LO amplitude, in order to achieve the exact RG in- 
variance. The primordial and distorted-wave counterterms are the short-range operators necessary 
to absorb the divergences of multiple pion exchanges sandwiched between free and LO interacting 
states, respectively. 

We have argued that RG invariance provides two mechanisms to enhance, relative to WPC, 
the short-range forces in ^Sq. (i) As the residual counterterms for the LO amplitude, they scale 



44( 1 : for a generic ^Sq counterterm 



similarly to the contact interactions of the pionless theory 
C2n with 2n derivatives, 

Att 1 

Cln • (35) 

(ii) As the distorted-wave counterterms for multiple pion exchanges, they are enhanced by 
©(M^j/Mj^), but only starting from the four-derivative term, 

Cfn ~ , \„ o , (36) 

Since the enhancement due to the residual counterterms dominates, we power count ^S'o contact 
interaction according to ([35|) . as if the theory were the pionless one. In particular, the 0{Q^) 
counterterms have the following structure: 

(i5o|yf 1^50) = + ^ (/^ + ^(^y V + ^P'V [p" +P') ■ (37) 

Because of the uncertainties associated with Cj in TPEl, we refrain from actually calculating 0{Q^) 
of ^5*0. Summarized in Table [I] is our power counting for the two-nucleon sector in both singlet and 
triplet channels for S and P waves. 

We scrutinize WPC with a more stringent interpretation of RG invariance: not only should 
the cutoff dependence become vanishingly small for A > Mhi, but it must vanish sufficiently fast 
so that the accuracy claimed by the power counting is consistent with the cutoff error. This 
leads to a crucial conclusion in our analysis that, contrary to WPC, 0{Q) of the EFT expansion 
does not vanish. Instead, 0{Q) is made of one insertion of the two-derivative ^Sq counterterm: 
D/2{p'^+p^). Although we are not the first to propose this, our argument, that the cutoff error 
of the LO amplitude is one order lower than TPEO and has to be corrected by the D term alone, 
provides some new insights. For instance, unlike Ref. Q] (also discussed later in Ref. 45|), our 
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0{1) 


OPE, 

V V 


„ , (Csp.p'po' 


) 


0{Q) 






TPEO, SiSoP'V, 
D^p„p'p{v'^ +P^), P'P 


ssDP'' y' 

^Dsp^{p'^+p^) EpT^p^] 

y EPFP'^ ) 

, Cap^p'p 




0{Q') 


TPEl, Fis,p'^pHp'^+P^) 



TABLE I: Power counting for pion exchanges, S and P-wave countcrtcrms up to 0{Q^). p {p') is the 
magnitude of the ccnter-of-mass incoming (outgoing) momentum. The two-by-two matrices are for the 
coupled channels. 



rationale is a priori and does not rely on the numerical value of -D in a particular renormalization 
scheme. 

A full, nonperturbative RG analysis, with OPE as the only long-range force, of the counterterms 
was attempted in Refs. [ll, [2^) which also concluded that the D counterterm is more important 
than TPEO. Although the nonperturbative RG analysis appears to be free of any guesswork for 
obtaining power counting, the robustness of the conclusions of Refs. [ill . [20I is obscured by the 
assumptions made therein to derive and solve the RG equation. On the other hand, our approach 
can be viewed as the explicit, order-by-order examination of an ansatz — the proposed power 
counting — to the RG equation. If RG invariance can be shown to hold at all orders, which 
we could not rigorously achieve though, we cannot think of any reason why the proposed power 
counting could not be one of the solutions to the RG equation. 

It is instructive to compare the power counting of with that of the attractive triplet channels. 
A non- vanishing 0{Q) arising in ^5o but not in the triplet channels has everything to do with the 
fact that OPE is regular (1/r) in ^S'o but singular (l/r^) in the triplet channels. It is interesting that 
the singular attraction of OPE costs a few more LO counterterms in the attractive triplet channels 
(e.g., ^Pq and — '^^2) but in the mean time it avoids the pionless theory-like proliferation of 
subleading counterterms. 

The distorted-wave enhancement to the singlet-channel short-range forces occurs in only S wave 
(^5o), and it affects the power counting to a lesser extent than that of the residual counterterms. 
In contrast, the distorted- wave enhancement in the attractive triplet channels takes place in higher 
partial waves (^Pq, ^P2 — ^F2, etc.) but not in S wave, and it plays more important role in power 
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counting than the residual counterterms. 
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Appendix A: LO wave function 

With the regularized LO potential, the S'-wave radial wave function is well defined at the origin 
and can be written as 

Mr) = MO)Mr) , (Al) 

where (pki^) is the regular solution in the sense (pkif) joikr) as r — )• 0, with jo{p) the zeroth 
spherical Bessel function. 

With regularization IV (r) — the Fourier transform of Vyiq) — becomes relatively flat on the 
inside while it resumes the Yukawa form on the outside. The LO contact potential, Vg'^ , is smeared 
inside and vanishes outside. This means that the inside wave function is largely decided by C^"-* (A) 
and A, whereas the outside part is dominated by a combination of the irregular (T-Lk{r)) and regular 
[Jkir)) solutions to the Yukawa potential, 

</.fc(r) = AA(C(°), A) [?^fc(r) + e(C(0), A)Jfc(r)] , r > A'^ . (A2) 

T-Lkir) and Jk{r) have the following small kr expansions: 

n=0 

-Hkir) = — V 5„ (K^r) {krf^ - 2Jk{r) In (K^r) , 

K — T ^ ^ 



(A3) 



n=0 
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where (,n{x) and Sn{x) are analytic functions around x = and can be further expanded to obtain 
the expansions of Jk{f) and ^fc(r) in powers of r. 

On the other hand, the three-dimensional (in-state) wave function is related to the LO T-matrix 

by 



A 



V, (f) = e^fc-+/ -^e-'- : • (A4) 



Therefore, V'fc(O) is given by 



Using Eqs. Q, ([JO]) and ([U]), we obtain 

A/" and in Eq. ()A2p are functions of C'^'')(A) and A because the inside and outside wave 
functions need to match near r ~ A~^ when A; = or any other low momentum where we decide 
to fit C^^\A). In the meantime, the behavior of ^kif) at r ^ 00 must be RG invariant in order 
for extracting scattering observables, 

= 7^ .^rmiy^ . + , r»A-i. (A7) 

[C('^>) —Ik 

Recalling that x(A;; A;)/[(C^'^^)~"'^ — Ik] is RG invariant and comparing the above wave function at 
r — 7> cx) to the form of an asymptotic state, we find that Af/C^^^ and 9 are both RG invariant. 

Appendix B: Useful integrals for subleading T-matrix 

We briefly describe the computations of the integrals stemming from insertions of subleading 
counterterms into the LO T-matrix, such as Eqs. (j2ip . (j26p . (j29p and ()30p . Note that the rather 
formal results shown here are not entirely independent of the regulator and we do not use the 
expressions derived here in the numerical calculations. However, they have enough mathematical 
rigor to illustrate how renormalization is carried out. 

When evaluating (1 + T(0)G)uJ^^(GT(°) + 1) 1^, we need to know the following integral 

^ ^ (2vr)3^ E-J^ + ie- ^^^^ 
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By differentiating with respect to x on both sides of Eq. ()A4p and letting x = in the end, one 
finds 



(B2) 



On the other hand, from the Schrodinger equation for the LO wave function, 

i-^V^Arl^) = W - mNV^^\x) \ ^f:{x) (B3) 
where V^^^ {x ) is the LO potential represented in coordinate space. With the above equations and 



Eq. (|A6p . one can eventually arrive at Eq. (j2ip . 

Successive differentiation on both sides of Eq. (|A4p and letting x = in the end relates integrals 
of the following form, 



(B4) 



to higher derivatives of the wave function at r = 0, which can in turn be related by successive 
differentiation on Eq. (|B3p to V^^^ {x ) and its derivatives at r = 0. 

Integrals involving the LO interacting Green function are encountered in computing Eq. (j26p : 



d^h dH2 



(27r)3 (27r) 



with 



4(/2,^i) = (2vr)= 



3 --l '2 'Jk{h, h) , 



+ 



E-^ + ie ' (£;--^ + ie)(£;--^ + ie) 
The generating function for these integrals is the coordinate space version of : 

dHi d% 



Its second derivative is given by 



(2^)3 (27r)3 



4(r2,ri)e*('2-^2-'i-^i). 



(B5) 



(B6) 



(B7) 



i-^^\2,i}'^k{x2.xi) = -mN5f{S2-xi)+ -mNV^''\x{2,i}) «4(x2,xi), (B8) 



where 



5f\x) 



dH 
(2^ 



Al-x 



(B9) 



Again, successive differentiation with respect to xi and/or X2 on both sides of Eq. (jBSP and letting 
2^1,2 = leads to the integrals (jBSp with larger m and/or n. 
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